Superradiant quantum phase transition in a circuit QED system: 
a revisit from a fully microscopic point of view 
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In order to examine whether or not the quantum phase transition of Dicke type exists in realistic 
systems, we revisit the model setup of the superconducting circuit QED from a microscopic many- 
body perspective based on the BCS theory with pseudo-spin presentation. By deriving the Dicke 
model with the correct charging terms from the minimum coupling principle, it is shown that the 
circuit QED system can exhibit superradiant quantum phase transition in the limit A'^ oo. The 
critical point could be reached at easiness by adjusting the extra parameters, the ratio of Josephson 
capacitance Cj to gate capacitance Cg , as well as the conventional one, the ratio of Josephson energy 
Ej to charging energy Ec- 
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Introduction. — In the conventional Dicke model 
which ignores two-photon interaction term A^, the su- 
perradiant quantum phase transition (QPT) can 
happen when the atom-radiation field coupling g is strong 
enough. However, for the realistic systems with the min- 
imum coupling from [/(l)-gauge theory, the Thomas- 
Reiche-Kuhn (TRK) sum rule means that the term 
cannot be neglected when the increases of A^ term fol- 
lows the increase of g 0| ■ The A^ term shifting the 
effective frequencies hence prevents the considered sys- 
tem from reaching the critical point, so that no superra- 
diant QPT happens in the natural atom systems. This 
fact was stated as a no-go theorem for the realistic cavity 
QED systems. 

On the contrary, it was ad hoc pointed out that this 
superradiant QPT could be realized in some artificial 
system[7(]. Later on, P. Nataf and C. Ciuti showed 
that the circuit QED system consisting of a collection 
of Josephson atoms capacitively coupled to a transmis- 
sion line resonator (TLR) is capable for such kind of QPT 
0. Viehmann et al questioned this judgment based on 
an overall microscopic model. They argued [9-11] that 
the phenomenological Hamiltonian used in Ref . t8j] cannot 
adequately describe the superradiant QPT of the circuit 
QED system with large atom numbers; if all the degrees 
of freedom are considered properly, the no-go theorem of 
superradiant QPT based on the TKR sum rule still works 
for this artificial system. In this sense, they excluded the 
existence of the superradiant QPT in the circuit QED 
system. 

Viehmann et al claimed that a fully microscopic ap- 
proach was utilized by themselves, but it seems diffi- 
cult to straightforwardly deduce the phenomenological 
Hamiltonian used in Ref. [8| from their overall microscopic 
model. In this paper, we try to carry out this necessary 
task to deduce it from a microscopic model with the min- 
imum coupling form U (l)-gauge theory. To this end, we 
provide a description of the circuit QED system using 



Figure 1: (color online) The schematic of the superconducting 
circuit for superradiant QPT. An array of Josephson junctions 
[with zoomed picture in (a)] is embedded in a superconducting 
TLR and couples to the resonator mode capacitively (b). The 
equivalent circuit for a single Josephson junction is depicted 
in (a), where the quantized resonator mode Vq is included in 
the gate voltage Vg. 



the pseudo-spin representation [1^ of the BCS theory, 
which explicitly displays the superconducting character- 
istics of the Josephson atoms. Our microscopic approach 
correctly gives the additional quadratic quantum voltage 
term, which is usually ignored in current references, e.g., 
[l^l- Applying this result to the low excited ensemble of 
artificial atoms, we conclude that the no-go theorem in 
the cavity QED system could not rule out the superra- 
diant QPT in the circuit QED system for some experi- 
mentally accessible parameters. 

Microscopic modeling of superconducting circuit 
QED. — As a key element in the circuit QED system, 
as illustrated in FiglU superconducting Josephson 
junction consists of a thin insulating barrier sandwiched 



2 



between two superconductors. Microscopically, we use 
the collective pseudo-spin operators 



2 H 



i) (1) 



to describe the charging (tunneling) process of Cooper 
pairs in the junction [T^. [l5| . Here, cjj.^(ck,Q) is the 
electron creation (annihilation) operator for the super- 
conductor on the a-hand side with a = L, R. The in- 
dex k (— k) denotes the momentum of the electron with 
spin up (down), and it is summed over the energy shell 
[ep — fy-^D,^F + fy^o], which is around the Fermi energy 
Ei? up to the Debye frequency ub ■ The normalization fac- 
tor Af equals to half of the number of momentum states 
within this energy shell. Obviously, 5^ „ is the operator 
counting the number of the Cooper pairs in excess of the 
electroneutrality of the superconductor on the a side. 

To model the tunneling process as the Josephson ef- 
fect [l2| , the single electron tunneling Hamiltonian Ht — 

-^X^k q ("^k fl'^q-i h.c.^ is re-expressed in terms of the 
collective pseudo-spin operators as 



Ht = -2hM^T {S+ + S- 



(2) 



where we have ignored the single electron tunneling terms 
for the system in superconducting phase. Here, the op- 
erator S~ = SlS]^ {S'^ = {S~y) denotes that a Cooper 
pair tunnels from the left (right) superconductor to the 
right (left), where the operator Sz — {Sz,l — Sz,r) /2 
is defined and the commutation relations [S'^jS'^] = 
TS± and [S+,S-] « 2Sz/JV^ are fulfilled. Note that 
Sz,L + Sz.B. ~ since the Josephson junction is elec- 
troneutral, Sz also represents the number of the excess 
Cooper pairs on the left bulk of superconductor. In the 
case of small number of the excess tunneling Cooper 
pairs, i.e., Af ^ (Sz), S'^S~ is central to the alge- 
bra generated by and Sz because of [S'+,S'^] 
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Then the polar decomposition S^ = exp (±z6') 
defines the macroscopic phase operator 6, which obeys 
[0, 5*2] — i. Using the phase operator, we can rewrite 
the tunneling Hamiltonian in the conventional fashion 
Ht = —Ej cos 9, where Ej = <I>o//27r is the Josephson 
energy, $0 = ^/2e is the quantized flux, and / = 8eN^T 
is the maximum tunneling current. 

According to the reference , the microscopic mean- 
ing oi 9 — (j)L — 4>R could be explained as the differ- 
ence between the order parameters of the right- and 
left-hand superconductors in the BCS ground states 

\'t>a) = Ilk "k + t'kexp(i0a)cLk,a4,a l^), whcrc 0q is 
the common phase of superconductor. Then the tun- 
neling current (,/) ~ /sin((^L — (pn) is obtained by the 
average of J = — 2e [A''^, 77^] over the product state 



\<t>) = I'Pl) "Si \4'r)7 where Nl is the number of the elec- 
trons on the left-hand superconductor which equals to 
2Sz plus a constant. 

Next we model the charging process for a simple 
Josephson device, which is a superconducting island [or 
Cooper pair box (CPB)] connected to a gate capacitor 
Cg and a bulk of superconducting electrode through a 
thin junction with capacitance Cj . The geometry of the 
superconducting circuit is shown in Fig[TJa), where the 
voltage Vc is applied to the gate capacitor by a classical 
source. This device is also coupled to a quantized electro- 
magnetic field provided by a superconducting TLR in a 
coplanar- waveguide geometry, which gives an additional 
quantum voltage Vq = Vq (^a + a^) where (a) is the cor- 
responding creation (annihilation) operator for the single 
mode of the TLR with eigenfrequency ujr- 

With the charges Q distributed on the island, the elec- 
trostatic potential V of the junction is determined by 



CjV-Cg{Vg-V)=' 



(3) 



with the total gate voltage Vg = Vc + Vq. Initially, we 
assume no excess electron exists, i.e., Q — 0, hence the 
potential of the CPB Vq = CgVg/ {Cg + Cj) is formally 
quantized as an operator. The total electrostatic energy 
U = Cg (Vo - Vgf /2 + CjVi/2 for both the gate and 
the Josephson capacitors connected to the electroneutral 
CPB is calculated as 



U^AEcnl^, 



(4) 



where Ec — e'^/2Cs is the charging energy for a single 
electron. Cs = Cj + Cg is the total capacitance and 
ng = CgVg/2e. 

After I excess electrons are added in the CPB, the to- 
tal energy is the electrostatic energy U plus the work W 
done to tunneling Cooper pairs. W is actually the work 
cost by the excess electrons to cross the barrier, which 
is actually supplied by the voltage source. The corre- 
sponding potential Vi is calculated by substituting the 
excess charges Q = Qi = — e^' ^Q rij into Eq.Q, where 

rij = cl^Cj is the single electron number operator. As we 
concern the charge accumulation process, the momentum 
states indexes k of the electron operators are of no im- 
portance. Instead, we assign to each electron operator a 
subscript j indicating the order of accumulating on the 
island. 

According to classical electrodynamics, to add one 
more electron on the island with I excess electrons al- 
ready on it, the work is calculated according to the for- 
mula Wi = —eni^iVi. This formular seems phenomeno- 
logical, but now we can derive it from the minimum cou- 
pling principle based on U (l)-gauge theory with a single 
particle Hamiltonian H^ ~ [p ~ eA (x)] /2m — e(j) (x). It 
describes an electron moving in the vector potential A (x) 
and scalar potential (j) (x) . Here we use the coulomb 
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rauge V • A fx) 



and the dipole approximation with 



A (x) w A (xo) = Ao and (x) sa V0 (x) ^=^0, which 
is consistent with the prerequisite of the discussion about 
the superradiant phenomenon in this paper. It leads to 
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(5) 



where Eo = - [dA (x) /dt + Wcj) (x)] 

In second quantization, the field operators ip (x) = 
X^k ''/'k (x) Ck is used with (x) approximately being 
the plane wave. Then the energy cost of a single elec- 
tron crossing the capacitor from one electrodes to the 
other one at d apart is calculated as 



(x) iJe (x) - (x + d) i/e^ (x + d) = erW, 



where V = d ■ and n = X^k'^k'^k- The momentum 
term and quadratic vector potential term Aq are both 
canceled out, and the remaining term —enV verifies our 
phenomenological formular of Wi . 

When TV excess electrons are added in the CPB, the 



total work W — J2fLo^ obtained as 

W = 4EcSl - SEcngS, - 2EcS, 



(6) 



where we use the fact Sz = '^j/^- "^^^ linear term 

2EcSz can be neglected because it will merely shift Ug 
by 1/4, which can be adjusted by tuning the gate voltage 
without influence the further discussion. 

At last, the charging Hamiltonian He = U + W is 
explicitly written as 



He = AEc {Sz 



(7) 



We remark that the last term /irtg = AEgjC j /Cg — 1)?!^ 
was neglected in some current references |l8l . Il9| , since it 
is a constant for the classical voltage and not related with 
the charges in the CPB. However, in the case of the gate 
voltage Vg contains a quantized component, this term 
provides a nonzero quadratic voltage term, which is ev- 
idently crucial in determining whether the superradiant 
QPT exists. 

Superradiant QPT in the Dicke model based on cir- 
cuit QED. — Now we further consider the circuit QED 
system, as shown in Fig[T]Jb), with a TLR coupled to 
N small junctions, which are modeled as the artifi- 
cial atoms of two energy levels. The total Hamiltonian 
Heir — '^7=1 Hj + hw^a^a is defined by 



H, = AEciSz 



E.J cos ( 



(8) 



which correctly includes the tunneling part Eq.(l2]) and 
the charging part Eq.([7]). Obviously, Heir is very similar 
to the cavity QED system for the atoms interacting with 
cavity modes through minimum coupling. Generally, it 
is very hard in experiments to realize the strong atom- 
field coupling in the conventional cavity QED systems. 



but the strong coupling regime is feasible in the current 
experiments of the superconducting circuit QED. There- 
fore, the circuit QED system is more ideal to investigate 
the superradiant QPT. 

At the degenerate point the Josephson junction be- 
haves as a two-level system, and the total Hamiltonian 
reads 



Hr 



Er 



EcCgVg 



N 

E 

+hajra^ a + hD (a^ + a)'^ + hF (a^ + a) , (9) 



where the correct two-photon term is included with HE = 
fiNCgVg/2e and HD = N EcC jC gV;} / . ctJ = (0| + 
|0)^. (1| and (j| = \{)). (0| - |1)^. (1| are defined in terms 
of the two lowest eigenstates |0)^. and of Sz,j- Here, 
we neglect a constant N (CjjCg — 1) /4 as it is a pure c 
number . 

To study the superradiant phenomenon in this circuit 
QED system, it is necessary to explore the circumstance 
that the atom number TV is large and the total excitation 
number is low. In this case , the collective excitation 
operator 



6''' = ^='> le) 



El 
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(10) 



behaves as bosonic operator in the atomic quasi-spin 
wave [231, which is defined by the eigenstates |e)^- = 

* (|0)^. - |1)^.) /^/2 and \g) ^ = (|0)^. |1)^.) /V2 of the 
j-th CPB. Then, H^ir is rewritten as 



Hr 



hi^rf^ a -\- hbjjb^b — ihQ (c 



t 



(6t - b) 



hEi 



(11) 



Apparently, our circuit QED system is reduced into an 
equivalent system of two coupled harmonic oscillators 
(CHO) with frequencies ujj = Ej/h and ujr, and the 
coupling strength is 



2VNEcCgVg 



he 



(12) 



Generally, we consider a CHO system with canonical 
coordinates xi and X2-, eigenfrequencies coi and lu2, and 
masses mi and m2, respectively. If the coupling term 
were inappropriately chosen as —gxiX2, the eigenval- 
ues of the coupled system would be imaginary when the 
coupling strength g is strong enough, specifically, when 
g^ > mim2U}fui2- Somebody depicts this phenomenon as 
a kind of QPT, but the natural coupling in the conven- 
tional coupled CHO should be g {xi — X2)^ /2 = gx\/2 + 
gx\l2 — gx\X2, so the two quadratic coordinate terms 
renormalize the eigenfrequencies as Cji = y^ojf + g/rrti, 
i = 1,2. Therefore, there would not be QPT, since 
g^ < TOiTO2Wi£2'| is always valid [21]. This is the very 
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reason that the correct non-hnear term is particularly 
important in the discussion of QPT. 

However, in our present circuit QED system, there is 
no such intrinsic relation between the renormalized eigen- 
frequency and the effective coupling, thus it is possible 
to observe such kind of QPT phenomenon wherein. We 
would like to point out our model Hamiltonian contains 
an additional linear term HF (a^ + a) , which is intro- 
duced accompanying the correct two-photon term. It can 
be eliminated by displaced transformations a = a + rj and 
/? = 6 -f ^ [r] and ^ are c numbers) . Then we diagonalize 
Heir hi a conventional way (2^ and obtain two eigenfre- 
quencies as 



(13) 



where r2± = y/uj^ + ADuir ± Wj. It is obvious that is 
always positive, while the can be negative when 



1 — K7 > 



UJrOJj 



(14) 



Here, we define two dimensionless parameters, the ratio 
of Josephson energy to charging energy k = Ej / AEc and 
the ratio of Josephson capacitance_to gate capacitance 
7 = Cj/Cg. It was proven in Ref.[3| that a;_ = is the 
critical point of the superradiant QPT and the superra- 
diant phase lies in the region that eigenfrequency uj- of 
the system is imaginary, hence Eq. (fT4l) is actually the 
condition for the appearance of superradiant phase. As 
the coupling fl increases with Vn, the right hand side of 
Eq. ()14p approaches to positive infinitesimal in the limit 
N ~^ 00. Thus the occurrence of the superradiant QPT 
depends on the condition K7 < 1. The corresponding 
critical point of f2 is 



1 /W^ 

2 V 1 - K7 



(15) 



According to the above arguments, the superradiant 
QPT indeed can occur in principle, but we need to ex- 
amine this conclusion for the realistic systems. In order 
to achieve a good charge qubit with small fluctuation of 
the Cooper pair number and low classical noise, it is usu- 
ally chosen k <C 1 and 7 ^ 1 in experiments. These two 
key factors, k and 7, compete in determining whether 
the superradiant QPT can take place. If the Josephson 
capacitance is large enough, such as 7 « 10 [S^l, and the 
ratio K is set around 0.4 [2J|, the condition Eq.(fT4| is 
violated thus the superradiant QPT can not happen. In 
contrast, we can also choose k ^ 1 and 7 « 1 as in the 
experiment |25j . which clearly allows the superradiant 
QPT. 

Remarks and conclusion. — Pedantically, we need to 
understand why the Dicke-type superradiant QPT is al- 
lowed in the circuit QED system, while it is forbidden in 
the cavity QED system, since these two systems possesses 



natural atom 


artincial atom 


H„a = pV2m + U (x) 


Har = 4Ec (5, - 1/2)" - Ej cose 


X 


e 


P 


- 1/2 


[x,p] = ih 


[d, - 1/2] = i 



Table I: The correspondence between the natural atom and 
the Josephson-type artificial atom. 



very similar Hamiltonians with the correspondences be- 
tween canonical variables as listed in the table HI This 
analogy apparently implies the superradiant QPT can 
happen nether in the circuit QED system nor the con- 
ventional cavity QED system. However, this argument 
obviously contradicts with the conclusion made above, 
as well with the analysis by Nataf et al [IC^ . 

To solve this puzzle, we would like to consider whether 
or not there exists the correspondence between the basis 
vectors used for defining the collective operators in two 
systems. In cavity QED system, we use the two lowest 
eigenstates |e) and \g) of the total Hamiltonian Hna of a 
natural atom to define a qubit subspace. Evidently, they 
are not the eigenstates of the momentum operator p due 
to the existence of trapping potential. In the artificial 
atoms, however, though the electron pair number oper- 
ator Sz corresponds to momentum p, the two discrete 
eigenstates |0) and |1) of Sz does not correspond to |e) 
and \g) respectively. It follows this observation that the 
collective operators of the natural atom ensembles and 
artificial atoms are of different types, and describe differ- 
ent types of quasi-excitations. Thus it is not surprising 
that the circuit QED system exhibits superradiant QPT 
while the cavity QED system does not. 

In summary, we have theoretically explored the super- 
radiant QPT in the circuit QED system, where A'' Cooper 
pair boxes behaves as artificial atoms coupled to a sin- 
gle resonator mode. With the microscopic Hamiltonian 
based on the pseudo-spin representation of the BCS the- 
ory and the minimum coupling principle, we deduce the 
correct quadratic term Vg of the gate voltage from a fully 
quantum perspective. Then we showed that the circuit 
QED system is capable for the superradiant QPT, and 
the critical point is determined by more parameters, the 
ratios k and 7. The QPT is more feasibly to be real- 
ized when these two ratios are small. We also explained 
the cavity and circuit QED systems show different collec- 
tive behaviors is due to the superradiant phenomenons in 
these two systems are based on different types of quasi- 
excitons. 
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